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We show that the recently introduced 0-1 test can successfully distinguish between strange 
nonchaotic attractors(SNAs) and periodic/quasiperiodic/chaotic attractors, by suitably choosing the 
arbitrary parameter associated with the translation variables in terms of the golden mean number 
which avoids resonance with the quasiperiodic force. We further characterize the transition from 
quasiperiodic to chaotic motion via SNAs in terms of the 0-1 test. We demonstrate that the test helps 
to detect different dynamical transitions to SNAs from quasiperiodic attractor or the transitions from 
SNAs to chaos. We illustrate the performance of the 0-1 test in detecting transitions to SNAs in 
quasiperiodically forced logistic map, cubic map, and Duffing oscillator. 



PACS numbers: 05.45.-a,05.45.Tp, 05.45. Pq 

A strange nonchaotic attractor(SNA) is consid- 
ered as a complicated structure in phase space. 
Such a complex structure is a typical property 
usually associated with a chaotic attractor. How- 
ever, SNAs are nonchaotic in a dynamical sense 
because they do not exhibit sensitive dependence 
on initial conditions (as evidenced by negative 
Lyapunov exponents). These attractors are ubiq- 
uitous in different quasiperiodically driven non- 
linear systems. While the existence of SNAs has 
been firmly established, a question of intense in- 
terest is how do these attractors morphologically 
differ from other attractors and how they can be 
characterized and quantified. Two specific char- 
acterizations, namely finite time Lyapunov expo- 
nents and recurrence plots were proposed ear- 
lier in the literature. However, in the present 
work we show that the recently introduced mea- 
sure of 0-1 test clearly distinguishes quasiperiodic 
motion, SNAs and chaotic attractors and also al- 
lows one to detect different dynamical transitions 
from quasiperiodic motion to SNAs and SNAs to 
chaotic attractors. 



I. INTRODUCTION 

A central problem in the study of determin- 
istic dynamical systems is to identify different types 
of asymptotic behavior of a given system and to un- 
derstand how the behavior changes as the system pa- 
rameter changes [![. The asymptotic behavior can be, 
for instance, a periodic oscillation, a quasiperiodic mo- 
tion, a SNA, or a random or a chaotic motion 0, Q. 
There has been considerable work in the past addressing 
how dynamical systems develop chaos from periodic or 
quasiperiodic motions and how they can be characterized 



by various diagnostic tools such as Lyapunov exponents, 
fractal dimemsion and recurrence analysis. Some among 
these like the Lyapunov exponents and fractal dimension 
also require substantial amounts of data, free of noise, 
in order to perform well. Recently, a much simpler test 
for the presence of deterministic chaos was proposed by 
Gottwald and Melbourne 0-0] • Their 0-1 test for chaos 
takes as input a time series of measurements, and returns 
a single scalar value either for periodic attractors or 1 
for chaotic attractors. The test does not require phase 
space recontruction. The dimension of the dynamical 
system and the form of the underlying equations are in 
general irrelevant, and the input is just the time series 
data. The test returns a unit value in the presence of 
deterministic chaos and zero otherwise. Recently, many 
works have been carried out in this connection, includ- 
ing reliability of the 0-1 test Q, characterization of noisy 
symbolic time series Q and efficient time series detection 
in FPU lattices @. However all these studies focused 
only on distinguishing between regular and chaotic mo- 
tions. In the present work we investigate and characterize 
the performance of the 0-1 test when applied to quasiperi- 
odically forced dynamical systems exhibiting SNAs. 

SNAs which are commonly found in quasiperiod- 
ically forced systems were first described by Grebogi 
et al [10] . These attractors are geometrically strange 
as they are properly described by a fractal dimension 
but the largest nontrivial Lyapunov exponent is nega- 
tive, implying nonchaotic dynamics. Briefly speaking, a 
SNA has a complicated (fractal) structure but does not 
have a positive Lyapunov exponent. Trajectories sepa- 
rate at rates that are slower than an exponential diver- 
gence filj |. SNAs have been reported in many physi- 
cal s yste ms such as the quasiperiodically forced pendu- 
lum [12| , biological oscillators 1131 . and quantum parti- 
cles in quasiperiodic potentials [1J], which are also re- 
lated to the Anderson localization in the Schrodinger 
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equation [15|. Also, these exotic attractors were con- 
firmed by an experiment consisting of a quasiperiodi- 
cally forced, buckled, magnetoelastic ribbon [l6j], in ana- 
log simulations of a multistable potential [17| , neon glow 
discharge experiment [Hj], chemical oscillations and 
various electronic circuits [13, Hlj]. SNAs can be quan- 
titatively characterized by a variety of methods, includ- 
ing the estimation of Lyapunov exponents and fractal 
dimension [IJ [22|, [23| and spectral properties 0, [Hj]. 
The geometric strangeness of the attractor can be mea- 
sured through indices such as the phase-sensitivity ex- 
ponent [26H28j . while the chaoticity properties can be 
studied by examining the finite-time Lyapunov expo- 
nents [13, The transition from quasiperiodic to 
chaotic motion via SNAs can be studied by both the 
variance of finite time Lyapunov exponents and the recur- 
rence plots [30|, HH . The quasiperiodic (torus) attractor 
and SNAs are characterized by the Lyapunov exponents 
with zero and negative values, respectively. However, for 
quasiperiodically forced dynamical systems, the values of 
the Lyapunov exponents are negative for both quasiperi- 
odic attractors(tori) and SNAs. In this case the Lya- 
punov exponents fail to detect transitions from quasiperi- 
odic dynamics to SNAs. Moreover, analytical methods 
such as the approximation of quasiperiodic forcing by a 
series of rational approximations, or the analysis of the 
phase sensitivity of the dynamics demand a knowledge 
of the underlying map or differential equations which is 
often not available. However in the present work for the 
first time it is reported that the 0-1 test clearly distin- 
guishes tori, SNAs and chaos, from the time series alone. 
In particular this test indicates a value '0' for torus, a 
value in between '0' and T' for SNA, and for chaos a 
value that tends to T'. 

The mechanisms and routes for the birth of SNAs, 
their characterization in different dynamical systems, and 
tools for distinguishing different dynamical transitions 
to SNAs were reported in [32Tj4lj ]. Among the differ- 
ent routes to SNAs, Heagy-Hammel, fractalization and 
intermittency transitions are ubiquitous in quasiperiodi- 
cally driven nonlinear systems. In the present paper, we 
employ the 0-1 test to characterize different dynamical 
transitions to SNAs, in particular the Heagy-Hammel, 
fractalization and intermittent transitions. Specifically, 
the 0-1 test characterization for transition to SNAs is 
examined in quasiperiodically driven logistic map, cubic 
map and parametrically driven Duffing oscillator to val- 
idate our conclusions. 

The organization of the paper is as follows. In 
Sec. II we briefly review the 0-1 test. In Sec. Ill we re- 
port the applicability of the 0-1 test to quasiperiodically 
driven systems. Sec. IV examines the characterization us- 
ing the 0-1 test for different routes of creation of SNAs, 
while in Sec.V we discuss the transition from SNAs to 
chaos. The final section gives our conclusions. In Ap- 
pendix A we briefly describe our procedure to identify 
an optimal value for the parameter c required for the 0- 
1 test. In Appendix B we briefly analyze the behaviors 
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FIG. 1: (Color online)Plots of (i) translation variable p versus 
q (ii) mean square displacement versus n for the logistic map 
(Eq.(ll)) corresponding to (a) periodic dynamics for e = 0.0, 
a = 3.5 and (b) chaotic dynamics for e = 0.0, a — 4.00. Here, 
the total length of the time series N = 5 x 10 4 with n = N/10. 
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FIG. 2: (Color online) Application of 0-1 test for the force-free 
case of (a) logistic map (Eq.(ll)) and (b) cubic map (Eq.(12)): 
(i) Bifurcation diagram, (ii) Asymptotic growth r&te(K) from 
0-1 test as a function of system parameter. Here N = 5 x 10 4 . 



of the translation variables p and q, and the asymptotic 
growth rate if as a function of the length of the time 
series N. In Appendix C we present the nature of at- 
tractors in the quasiperiodically forced cubic map using 
the 0-1 test. Finally, in Appendix D we analyze the dis- 
tribution of local asymptotic growth rate P(K, N) along 
with finite time Lyapunov exponents (FTLEs) for a sta- 
tistical evaluation of the results. 



II. UNDERSTANDING OF 0-1 TEST 

Consider an n-dimensional dynamical system X = 
F(X), X = (xi, x 2 , x n ) T and denote a solution of the 
underlying system as X(t). Let (f>(x) be an observable of 
the underlying system. The 0-1 test uses the observable 
to drive the dynamics on a well chosen Euclidean exten- 
sion, and then exploits a theorem from Nicol et al [42ll43| . 
which states that the dynamics on the group extension is 
bounded if the underlying dynamics is nonchaotic, but 
it behaves like a Brownian motion if the dynamics is 
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chaotic. The 0-1 test for chaos first appeared in [4{. It 
is designed to distinguish chaotic behavior from regular 
behavior in deterministic systems. The test results in a 
value 1 if the system is chaotic and a if the system 
is regular. Theoretical justification was given to the 0-1 
test in |6|. An implementation technique was also pro- 
posed for the 0-1 test in which has been very useful 
for the present work. Here we present a brief outline of 
the 0-1 test followed by a modified version of the 0-1 test 
illustrated in dSl. 



A. Translation variables and the 0-1 test for 
regular and chaotic motion 

Consider a time series of length TV. The 0-1 test begins 
with a computation of two variables p{n) and q(n) which 
are the so-called translation variables: 



p(n) 



^2<t>(j)cos(jc), 

3 = 1 



<l( n ) = ^2^U)sin(jc), 

3 = 1 



(1) 



(2) 



where n — 1,2, N. 

Here <f>(J) is an observable constructed from the time 
series. The value of c may be chosen randomly in the in- 
terval (0,27r) for nonquasiperiodic systems @. One may 
note that it is the behavior of these translation variables 
that determines whether the dynamics of the system is 
regular or chaotic. The translation variables also exhibit 
a Brownian motion as a function of n, when the dynamics 
is chaotic. For Brownian motion we expect the average 
°f [pU + n )-p(j)} an d [q(j + n)-q(J)] to grow like y/n for 
j = 1,2,3,4... Therefore the square of [p(j + n)-p(j)] 
and [q(j + n)-q(J)] should asymptotically approach a lin- 
ear growth proportional to n for sufficiently large values. 
Hence, to find the behavior of the translation variables 
the 0-1 test calculates their mean square displacement 
(D(n)), 

1 N 

D(n) = Jim — V [p(j +n)- p(j)] 2 +[q{j + n) - q{j)) 2 

AT— >oo iV A — ' 

3 = 1 

However, in order to obtain a better convergence of mean 
square displacement, for numerical purpose, we have used 
a modified mean square displacement M(n) (which ex- 
hibits the same asymptotic growth rate as D{n)) of the 
form as suggested in Q 

M(n) = D(n)-V osc (c,n), (4) 
(1 — cos(nc)) 



where 



The theory behind the 0-1 test shows that if the dynamics 
of the system is regular, then (4) is bounded in time. 
However if the dynamics of the system is chaotic then (4) 
scales linearly in time. For numerical purpose, with the 
limit in (4), we must use n << N. In practice, the limit 
M{n) exists when n = N/10, where one expects M{n) 
to scale linearly. The dynamics is further characterized 
through the computation of the asymptotic growth rate 
K which can be calculated directly from M(n) by using 
a linear regression for the log- log plot of M(n) given by 
the defnition, 



K = lim 



logM (n) 
log(n) 



(5) 



It has been pointed out in [6j that the linear regression 
method leads to a distortion for small values of n. Hence, 
it has been recommended [6( to estimate the asymptotic 
growth rate K by using a correlation method which works 
better than the linear regression method. In the present 
work we employ the correlation method for the estima- 
tion of K. In this procedure, K is the correlation coeffi- 
cient of the vectors 



£= (l,2,....,n) T , 
and 

A = (Af(l),M(2),....,M(n)) T . 

Now one can define the correlation coefficient 

cov(Z, A) 



K = corr(^, A) 



6 -1,1 



(6a) 
(6b) 

(7) 



Here 



s/ i>ar(£)uar(A) 
1 " 

(£,A) - -X;«(j)-0(A(i)-A),«ar(0 = aw(f,0, 



3 = 1 



where 



1 ™ 



3 = 1 



(1 — cose) 



3 = 1 

Finally a value of K = indicates a non-chaotic data 
set while a value of K = 1 indicates a chaotic data set. 
We use the above analysis to characterize SNAs from 
quasiperiodic and chaotic attractors by properly choos- 
ing the value of 'c' in Eqs.(l)-(2) for quasiperiodic dy- 
namical systems. 

The functions pin) and q{n) in Eq.(l) can be viewed 
as a component of the solution to the skew product sys- 
tem, 

9{n + l) = 9(n)+c + <p(n), (8) 



p(n + l) =p(n) + (f>(n)cos(6(n)), (9) 



q(n + 1) = q(n) + <p(ri)sin(6(ri)), (10) 
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driven by the dynamics of the observable 4>{ri). Here 
(6, p, q) represent the coordinates on the Euclidean group 
rotation 9 and translation (p, q) in the plane. It was 
argued that an inspection of the dynamics of the (p,q) 
trajectories provides a quick and simple visual test of 
whether the underlying dynamics is regular or chaotic [H, 

a. 



B. Extending the 0-1 test to quasiperiodically 
forced systems exhibiting SNAs 

The above 0-1 test, when applied in the standard 
way to quasiperiodically forced systems, exhibiting 
SNAs, does not clearly distinguish SNAs from torus 
and chaotic attr actors. The 0-1 test returns either a 
value or 1 for the SNAs depending on the route by 
which they are created, as long as the parameter c in 
Eqs.((l)-(2)) is chosen as an arbitray number in the 
interval (0,27r). In fact, in a recent study Dawes and 
Freeland [44[ have attempted to distinguish between 
quasiperiodic attractors and SNAs with a modification 
of the 0-1 test by the addition of noise with the mean 
square displacement. In particular, by choosing the 
value of c to be equally spaced in the interval (0.1,tt/2), 
they showed that for a quasiperiodic attractor K takes a 
value zero and for a SNA it takes a value '1' with large 
deviations for certain values of c. However this method 
does not distinguish SNAs from chaotic attractors. In 
the present study, we have chosen the value of c in 
terms of the golden mean number like (s/Z + l)/2 or 

£^ , where u) = (\/5 — l)/2 to characterize SNAs in 
quasiperiodically forced dynamical systems, so as to 
avoid resonance with the translation variables. Details 
are given in Appendix A. With such a choice, we find 
that we can clearly distinguish torus, SNA and chaos as 
well as different dynamical transitions involving SNAs. 
These are discussed in detail in the following. 

In our numerical analysis the value of the asymptotic 
growth rate K is computed by taking the average of 
it over a number of segments in the given trajectory. 
That is, a given trajectory is divided into a number of 
segments each of length L and, in each of these, the 
corresponding asymptotic growth rate Ki (where V 
denotes the segment) is evaluated. Finally the average 
is taken to compute the value of K. Specifically, we 
consider a total time series of length iV = 3 x 10 5 for 
both quasiperiodically forced logistic map and cubic map 
while for quasiperiodically forced double well Duffing 
oscillator we take N — 5 x 10 5 (for continuous systems 
the data points are considered after sampling the time 
series data). In the total time series L = 10000 data 
points is considered as the length of each segment. 

We also note here that the qualitative nature of the 
attractors, the trajectories in the translation variables 
(p, q) space, and the mean square displacement essen- 



tially remain unchanged as a function of N (as discussed 
in Appendix B). So whenever these quantities are dis- 
played in the following, we use N = 5 x 10 4 data points, 
but for the calculation of K we use a much larger time 
series as mentioned above. 



III. APPLICABILITY OF 0-1 TEST IN 
QUASIPERIODICALLY DRIVEN SYSTEMS 

In this section we consider the applicability of the 
0-1 test for quasiperiodically forced dynamical systems 
to identify SNAs and distinguish them from periodic, 
quasiperiodic and chaotic attractors. We choose the 
following models of quasiperiodically forced dynamical 
systems for our present study; 

(i) Quasiperiodically forced logistic map: 

The dynamics is described by the coupled map [22j], 

x n+ i = a[l + ecos{2n9 n )]x n {\ - x n ), (11a) 
n +i =6 n + uj(modl), (lib) 

where u> = (vo— l)/2 is the irrational driving frequency, 
e represents the forcing amplitude, and a is a driving 
parameter. The different dynamical regimes can be 
characterized by using the rescaled parameter of the 
forcing amplitude e' = e/(4/a — 1) versus a. In this 
map the transition from torus to SNAs is characterized 
by the behavior of the largest Lyapunov exponent as 
well as the characteristic distribution of the finite time 
Lyapunov exponents. Further, SNAs are created at 
a saddle node bifurcation, when the dynamics shows 
type-I intermittency and other transitions. 

(ii) Quasiperiodically forced cubic map: 

The system is given by 

x n +i = Q + fcos(2ir9 n ) - Ax l + xf, (12a) 
n+1 = 9 n + ui(modl), (12b) 

where u> = — l)/2 is the irrational driving frequency. 
By varying the values of A and /, the creation of SNAs 
has been shown to occur by means of four distinct routes 
due to the truncation of torus doubling bifurcations. 
In particular, the formation of SNAs through Heagy- 
Hammel, fractalization, type-Ill intermittency and crises 
induced intermittency are described in [3l|. The SNAs 
created through different mechanisms are characterized 
in terms of the Lyapunov exponents, phase sensitivity 
and finite time Lyapunov exponents. 

(in) Damped two frequency parametrically 
driven double well Duffing oscillator: 

The damped two frequency parametrically driven double 
well Duffing oscillator is described by the nonlinear 
differential equation 
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FIG. 3: (Color online) For the quasiperiodically forced logistic map (a) Torus motion for e' = 0.3, a = 3.4874 (b) SNA for 
e = 0.3, a = 3.488 and (c) chaotic attractor for e = 0.5, a = 3.6: (i) Projection of the attractor; (ii) Dynamics of the translation 
components (p, q) in terms of 0-1 test (Eqs(l)-(2)); (iii) Mean square displacement as a function of n. Here N — 5 x 10 4 with 
n = N/10. 
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at 

which is a well suited model for buckled beam oscilla- 
tions. Here the driving frequency a is chosen to be the 
irrational number er = (v5 + l)/2. As the parameters 
R and H are varied one can locate different regions in 
the parameter space which exhibit different dynamical 
states. There are several mechanisms through which 
SNAs are formed in this system, two of which appear 
to be quasiperiodic analogies of intermittency transition 
in unforced systems. In addition to these, torus collision 
route to SNA as well as fractalization route have been 
identified in [23| . 



A. Dynamics in the absence of quasiperiodic force 

To start with we characterize periodic and chaotic 
attractors in Eqs.(ll) and (12) in the absence of 
quasiperiodic forcing. For this purpose the value of c in 
Eqs.(l) and (2) is randomly chosen in the interval (0,27r), 
and the observable is (f>U) = x j- First, the periodic and 
chaotic behaviors are investigated in the logistic map by 
using the new coordinates of translation variables p(n) 



and q(n). To analyze the behavior of the translation 
variables, let us fix the parameters as a — 3.5 (e = 0.0) 
in (11). Then the system admits regular dynamics. For 
a = 4.00 (e = 0.0) the system exhibits chaotic dynamics. 
For regular/chaotic dynamics the trajectories in the 
(p(n), q(n)) space indicate bounded/unbounded (Brow- 
nian) motion [see Fig. la(i), Fig. lb(i)], respectively. 
The bounded and unbounded behaviors in the {p(n), 
q(n)) space is further investigated in terms of the mean 
square displacement. Figs, la(ii) and b(ii) indicate 
mean square displacement versus n which is oscillatory 
for periodic dynamics, and scales linearly with time for 
chaotic dynamics. A similar behavior has also been 
confirmed for the cubic map with / = in Eq.(12). 

Next, we consider the full dynamics of the logistic 
map given by Eq.(ll) and cubic map (Eq.12) in the ab- 
sence of quasiperiodic forcing. While the logistic map 
exhibits period-doubling bifurcations as shown in Fig. 
2a(i), the cubic map exhibits period doubling route to 
chaos followed by period halving phenomenon as a func- 
tion of the system parameter Q which is shown in Fig. 
2b(i). In order to confirm the applicability of the 0-1 test, 
we evaluate the asymptotic growth rate K of the mean 
square displacement given by Eq.(7) in Figs. 2a(ii) and 
b(ii), respectively, for the logistic map and cubic map 
in the absence of quasiperiodic forcing. As expected the 
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FIG. 4: (a) Variation of the largest Lyapunov exponent (b) 
variation of the asymptotic growth rate (K) from 0-1 test as a 
function of a for fixed e — 0.05 indicating torus— »SNA— »chaos 
transitions in the quasiperiodically forced logistic map. The 
vertical arrows in this and the following figures indicate tran- 
sition points. Here, N = 3 x 10 5 with L — 10 4 as the length 
of each segment. 



0-1 test returns a value for periodic motion and 1 for 
chaotic motion. A similar analysis has been performed 
for the Duffing oscillator (13) with single periodic force, 
but we do not present the results here for brevity. 



B. Dynamics with quasiperiodic force 

For quasiperiodically forced discrete driven systems as 
well as flows, the dynamics can be chaotic, nonchaotic 
or strange nonchaotic depending on the values of the 
system parameters. Previous studies show that in 
quasiperiodically driven dynamical systems the re- 
gions of SNAs can be quite complicated and also the 
boundaries between chaotic attractors and SNAs are 
convoluted, with regions of torus attractors interspersed 
among the regions of SNAs [22|, [3j|. We now wish 
to characterize SNAs in these quasiperiodically forced 
systems by using the 0-1 test to distinguish SNAs 
with quasiperiodicity and chaos. Again, we take the 
time series data (observable) as </>(j) = Xj and the data 
always contains a component of the quasiperiodic forcing. 

In the present study, as noted earlier we fix 
a single value for c in Eqs.(l) and (2) in terms of 

an irrational number like (\/5 + l)/2 or , where 
w = (Vb— 1)/2 to characterize SNAs in quasiperiodically 
forced dynamical systems. We zero down on a most 
suitable value of c by working on the attractors in the 
space of the translational variables for a set of values of 
c in terms of the golden mean for typical cases of torus, 
SNAs and chaos in the quasiperiodically forced systems 



(11)-(13). We pick up that value of c which most closely 
replicates the expected nature of the attractors: torus 
like structure for quasiperiodic attractors, bounded 
motion for SNAs and unbounded motion for chaotic 
attractors. The details are given in Appendix A. By 



this procedure we find that a choice c = 



(2tt) 



where 



uj = (V5 — 1)/2, is the most suitable value for the present 
study. 



The different possible dynamical behaviors are 
characterized in the logistic map, by varying the values 
of e' = e/(4/a — 1) and a. For instance, for e = 0.3, 
a = 3.4874 the system admits torus like motion; for 
e' = 0.3, a = 3.488 it is a SNA, while for e' = 0.5, 
a = 3.6, it is a chaotic attractor [22|]. Figures [3] a(i),b(i), 
and c(i) clearly indicate the projection of the dynami- 
cal attractors in the (x, 0) space. Then, the dynam- 
ics of translation variables in the (p(n),q(n)) space in- 
dicates that the points are uniformly distributed over a 
circle for quasiperiodic (torus) regime (Fig.3a(ii), Brow- 
nian like behavior for chaotic regime (Fig.3c(ii)), and 
for SNAs the (p(n),q(n)) space shows that the points 
are uniformly distributed and bounded with minimal 
Brownian like behavior (Fig.3b(ii)). We characterize 
the behavior of the translation variables by using the 
mean square displacement(M(n)) as a function of n, and 
find that M(n) is oscillatory for quasiperiodic regimes 
(Fig.3a(iii)); grows linearly with time for chaotic attrac- 
tors (Fig.3c(iii)), and for SNAs M(n) grows linearly ini- 
tially with oscillatory motion, but tends to get bounded 
for large n (Fig.3b(iii)). The aforementioned analysis 
clearly shows that the 0-1 test allows one to distinguish 
SNAs from quasiperiodic regime and chaos. A similar 
analysis has also distinguished tori, SNAs and chaotic 
attractors in the quasiperiodically forced cubic map rep- 
resented by Eq.(12). For details, see Appendix C. 

We next study the characterization of different dy- 
namical attractors as a function of the parameter a in 
the quasiperiodically forced logistic map (Eq.ll) in the 
parameter regime a £ (3.33, 3.41) by using the 0-1 test. 
The logistic map (Eq.ll) has been already investigated 
and various transitions to SNAs have been identified us- 
ing tools such as Lyapunov exponents and their variance, 
finite time Lyapunov exponents and power spectral mea- 
sures [12, Hy. Here, we consider the above parameter 
regime investigated in Ref. [22j for our 0-1 test analysis 
also. To compare our results using the 0-1 test with the 
already reported results on torus — s-SNA— >■ chaos tran- 
sitions, here we present the nature of the largest Lya- 
punov exponent along with the asymptotic growth rate 
K using the 0-1 test as shown in Fig. 4. As noted above 
our present study of the 0-1 test using the value of c = 

£^ , where uj = (V5 — l)/2, indicates that the asymp- 
totic growth rate returns a value for torus regime, 1 for 
chaotic regime, while in the SNA regime the value of K 
varies in between and 1 non-monotonically. It may be 
also noted that there is an abrupt transition in the value 
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of K from a small value to a larger one (indicated by a 
vertical arrow below a = 3.37) during the transition from 
torus to SNA. Between the value of a, 3.35< a ^3.37, 
the torus gets wrinkled before transition to a SNA which 
is indicated by a nonzero value of K. One also notes 
that while the Lyapunov exponent (A) as a function of a 
(Fig. 4(a)) does not clearly distinguish between torus and 
SNA, where it takes negative value in both the regimes, 
the asymptotic growth rate K clearly distinguishes the 
torus and SNA regimes, as well as the chaotic regime. 

Similar analysis has been carried out for the 
quasiperiodically forced cubic map and Duffing oscilla- 
tor to confirm that with the choice of c = - — — , where 
uj = (\/5 — l)/2, as an irrational number, the asymptotic 
growth rate K varies between and 1. 

IV. DIFFERENT ROUTES TO SNAS 

One of our main concerns in this work is whether the 0- 
1 test can distinguish different dynamical transition sce- 
narios involving SNAs through appropriate characteriza- 
tions. In this section we show clearly that the forma- 
tion of SNAs through different routes exhibits distinctly 
clear signatures in terms of changes in the growth rate of 
K. In particular we demonstrate this for the ubiquitous 
routes, namely Heagy-Hammel route on torus doubling 
to chaos via SNAs [241 ], fractalization of torus leading to 
SNAs [111 H3] , and the appearance of SNAs through type 
-I intermittent phenomenon (22[. For other routes also, 
one can identify similar signatures. We also note here 
that each one of the routes has specific significance as dis- 
cussed in the literature [30l - l33l | : (i) Heagy-Hammel (HH) 
route: Collision of doubled torus with its unstable parent 
leading to SNA before the onset of chaos, (ii) Fractaliza- 
tion route: Gratual fractalization of torus giving rise to 
SNA before the onset of chaos, (iii) Intermittency route: 
Analog of saddle node bifurcation leading to the birth of 
SNA before the onset of chaos. 



A. Heagy-Hammel Route 

The first route we consider is the Heagy-Hammel(HH) 
route. In this route, a SNA occurs when the period 
doubled torus collides with its unstable parent. That 
is, the period-2" torus gets wrinkled and upon collision 
with its unstable parent, a period-2 n_1 band SNA is 
formed. To exemplify the nature of this transition, 
let us consider the quasiperiodically forced logistic 
map (Eq. 11), where we fix the parameter e' = 0.3, 
while a is varied. In this case, the HH transition 
takes place at othh — 3.487793 [22j. When we examine 
the neighborhood of othh, the transition from torus to 
SNAs is characterized by the Lyapunov exponent(A) as 
shown in Fig. 5a(i). Particularly in the torus region, 
ot ^ othh, A varies smoothly whereas in the SNA region, 
ol cthh, the variation in A is irregular. To validate 



the applicability of the 0-1 test for SNAs, we consider 
the aforementioned parameter region and analyze the 
asymptotic growth rate K versus a. One can see that in 
Fig. 5a(ii) K becomes minimal, that is '0' or nearly '0', 
in the torus region. When a othh, the growth rate K 
is varying between 0.2 and 0.7 non-monotonically. The 
behavior of K versus a shows a transition from regular 
dynamical regime to a strange nonchaotic one with a 
threshold value othh, separating both regimes. 



Similarly, in the quasiperiodically forced cubic 
map [3l|| given by Eq. (12), we fix the system parameter 
at / = 0.7 and vary A in the range 1.8865 to 1.8875. 
Initially for A ~ 1.8865 the system has a period doubled 
torus. Then as A increases to the value 1.8868, the 
torus begins to wrinkle and approaches its unstable 
parent. After this the torus collides with the parent and 
gives birth to a fractal attractor (SNA), when A reaches 
the value 1.88697. For this value of A, the attractor 
possesses a geometrically strange and nonchaotic be- 
havior [311 ]. It is noticed that these dynamical changes 
are clearly distinguished in the 0-1 test as shown in 
Fig. 5b(ii), while the corresponding Lyapunov exponent 
structure is shown in Fig. 5b(i). For the torus region 
the value of K is close to '0', and in the SNA region 
it oscillates between 0.2 and 0.4. We can again note 
an abrupt transition in the value of K at the transition 
point A = 1.88697 when the torus gets transferred to a 
SNA. 



We conclude that the results of 0-1 test indicate in 
both the quasiperiodically driven logistic and cubic maps 
that the value of K is smooth in the torus region and, at 
the critical value, there is a sudden jump after which K 
oscillates with irregular behavior. 

The HH analysis has also been carried out for the 
parametrically driven Duffing oscillator (Eq. (13)) with 
the system parameter fixed as R = 0.3, while the value 
of H is varied between 0.181 to 0.186. Here one finds 
that the transition from a SNA to a torus takes place 
at H c = 0.184 [23j]. The Lyapunov exponent A does 
not show any distinct change in its behavior during this 
transition. Figure 6(a) is a plot of A as a function of H, 
in the neighborhood of H c , where the variation of A is 
small and smooth, but in the SNA phase the variation 
of A is large. On the other hand, the plot between the 
asymptotic growth rate K and the parameter H indicates 
as shown in Fig. 6(b) that the value of K sharply (but 
smoothly) decreases from 0.78 to 0.00 within a narrow 
range of the parameter H with respect to the transi- 
tion from SNA to the torus regime, again giving a clear 
distinction between regions of SNAs and quasiperiodic 
attractors. We may also note that the HH transition is 
quite narrow in the parameter range H as seen in Fig. 6. 
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exponent(A) (b) Behavior of the asymptotic growth rate (K) 
from 0-1 test. Here, N — 5 x 10 5 (after sampling the data 
points) and L = 10 4 . 



B. Fractalization Route 



transition from a torus to a SNA (through a wrinkled 
torus (WT)) which occurs at a F = 2.6526 The 
Lyapunov exponent A is smooth and linear with no 
particular signature observed during this transition. 
However, as shown in Fig[7^(ii) the results of the 0-1 
test indicate that the transition from tori to SNAs 
via WT is characterized by the asymptotic growth 
rate K which grows linearly from torus to WT, and be- 
comes irregular at the transition for SNA at qj? = 2.6526. 

To study the fractalization mechanism in the cubic 
map(Eq.(12)), we fix the parameter at / = 0.1, and vary 
A in the range 2.14 to 2.17. For A = 2.14 the system 
admits a torus like motion; For A — 2.16 the system 
exhibits quasiperiodic oscillations of the doubled torus 
and when it is increased to A — 2.165, the doubled torus 
becomes a WT, which persists up to A — 2.167, where 
the oscillatory structure(wrinkled) gradually approaches 
a fractal structure(SNA) [3l|. This region is examined 
by both the Lyapunov exponent(A) and the 0-1 test. 
Again the Lyapunov exponent A does not show any 
specific signature when the transition takes place [See 
Fig. 7b(i)]. But in the case of the 0-1 test as shown 
in Fig. 7b(ii) the value of asymptotic growth rate K is 
nearly in the torus region and is linearly increasing 
during transition from torus to WT, and exhibits an 
irregular behavior at the transition to SNAs. 



The fractalization route is a common and rapidly 
changing intriguing transition to SNAs as pointed out 
by Nishikawa and Kaneko [26[. Here a period n torus 
attractor gets wrinkled and eventually forms a n-band 
SNA. Now we examine the fractalization route in the 
quasiperiodically forced logistic map (Eq. (11)), where 
we fix e' = 1.00 and the value of a is varied from 2.64 
to 2.66. Here, the variation of A as a function of the 
control parameter a as shown in Fig. 7a(i) describes the 



The fractalization route occurs in the parametri- 
cally driven Duffing oscillator with R = 0.3, and H in 
the range H = (0.072, 0.073). The transition from a SNA 
to a torus occurs at Hp = 0.0729 [23j. We find that at 
this transition value of H the asymptotic growth rate K 
decreases to zero from the SNA region to the torus one 
[see Fig. 8(b)], whose Lyapunov exponent is as shown in 
Fig. 8 (a). Again during the transition there is an irregu- 
lar behavior in K during transition from SNAs to WT, 
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10 s (after sampling the data points) with L — 10 4 . 



while it is fairly smooth and linear (with small fluctua- 
tions) in K from WT to torus. One can conclude that 
the asymptotic growth rate K behaves differently in the 
fractalization route compared to the HH mechanism. 



C. Intermittency Route 

Finally we discuss the intermittency route for SNA, 
which differs from those formed through other mech- 
anisms. In [22j the mechanism for the route in the 
logistic map has been studied. It has been shown that 



in this route a SNA disappears and is replaced by a 
one-frequency torus through an analogy with the saddle 
node bifurcation. In the logistic map [Eq. (11)] the 
intermittent transition from a SNA to a torus occurs at 
a c = 3.4058088 with e' = 1.00 [22]. The intermittent 
transition from a SNA to a torus in the cubic map(Eq. 
(12)) is studied in[3l|. For the forcing amplitude / = 0.7, 
a SNA occurs at A = 1.801685; On increasing the value 
of A, an intermittent transition occurs at A c = 1.8017 
which indicates a torus motion. 

In Fig. 9 (a) we present the nature of the Lya- 
punov exponent in the neighborhood of the transition 
point in the system parameter both for the logistic 
map (Fig.9a(i)) and cubic map (Fig.9b(i)). At the 
intermittent transition, the Lyapunov exponent does not 
show any significant variation during the transition from 
SNAs to torus. On the other hand the 0-1 test detects 
intermittent SNAs and torus clearly in Figs. 9a(ii) and 
b(ii), where the asymptotic growth rate K fluctuates in 
the intermittent SNA region while it smoothly varies 
in the torus region. The value of K exhibits a sharp 
transition from a positive value (nearly 0.8) to zero 
during the transition from SNAs to torus. 

Further we consider the intermittent transition in 
the parametrically driven double well Duffing oscillator 
[Eq. (13)] with R = 0.35, while the value of H is varied. 
For H = 0.19080564, the system admits a chaotic attrac- 
tor; a SNA appears at H — 0.190833. Further when the 
value of H is increasing, an intermittent transition from 
a SNA to a torus occurs at H c = 0.19088564 [Hj]. At 
this transition point, the Lyapunov exponent [Fig. 10(a)] 
is of decreasing nature, where for both intermittent SNA 
and torus the value of A is negative. The 0-1 test clearly 
detects changes in the growth rate as shown in Fig. 10(b). 
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FIG. 11: (a) Variation of the Lyapunov exponent A and (b) 
variation of the asymptotic growth rate (K) from 0-1 test as 
a function of a indicating transition from SNAs to chaotic 
attractors in the quasiperiodically forced logistic map with 
e = 0.30. Here N = 3 x 10 5 and L = 10 4 . 



V. THE TRANSITIONS FROM SNAS TO 
CHAOTIC ATTRACTORS 

In the previous sections we have characterized different 
routes for the transition between tori and SNAs, which 
are qualitatively distinguished in terms of the asymptotic 
growth rate K by the 0-1 test. It has been shown that the 
value of K shows an abrupt change followed by an irreg- 
ular behavior in general after the birth of SNAs. In this 
section we mainly concentrate on important dynamical 
characteristics of transition from SNAs to chaotic attrac- 
tors or vice versa, usually characterized by the largest 



Lyapunov exponent which changes from negative to pos- 
itive values (22), IH, |36| . In general SNAs occur in nar- 
row parameter ranges, and finally get transferred into 
chaotic attractors. The transition is not usually accom- 
panied by any major change in the form or shape of the 
attractor. We now study this type of transition in the 
quasiperiodically forced logistic map [Eq.(ll)] with pa- 
rameter e' = 0.3. Negi etal [4l[ have found that the 
transition from SNAs to chaos occurs at a = 3.512. Fig. 
11(a) presents the Lyapunov exponent (A) for the tran- 
sition from SNAs to chaos which clearly shows that A 
changes from negative to positive values. In the case 
of the 0-1 test, the value of the asymptotic growth rate 
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(nearly 0.8) to zero during the transition from SNAs to 
torus. Finally, our analysis shows that for quasiperiodic 
attractors the value of K is zero, and if chaotic it is 1, 
whereas for SNAs it takes values in between '0' and '1'. 

Thus the tools of 0-1 test helps one to distinguish 
SNAs from other attractors in a clear manner. Further, 
we note that the advantages of 0-1 test are that it requires 
minimal time series data, does not require any equation 
of motion and is fairly easy to compute. This test has 
some disadvantages in that it depends upon a properly 
chosen value of V in Eq.(2). In this paper we have used 
'c' in terms of suitable irrational numbers to character- 
ize quasiperiodically forced dynamical systems. Finally, 
as it requires minimal time series, we believe that this 
test is more suitable for experimental situations in or- 
der to distinguish SNAs from other regular and irregular 
behaviors. 
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K oscillates between 0.5 and 0.8 in the SNA regime, and 
suddenly goes to T' in the chaotic regime [See Fig. 11(b)]. 
Similarly we have also carried out a study of the tran- 
sition from chaos to SNAs in the parametrically driven 
Duffing oscillator [Eq.(13)] with R = 0.47 while varying 
H. For H = 0.05, the attractor is chaotic. As H is in- 
creased to 0.0558, one obtains a SNA. The 0-1 test shows 
a specific signature of transition from chaos to SNAs as 
shown in Fig. 12(b) whose Lyapunov exponent is shown 
in Fig. 12(a). 

VI. CONCLUSION 

In this work we have studied different routes to SNAs 
and transitions from tori to SNAs and then to chaotic at- 
tractors in prototypical examples, namely the quasiperi- 
odically forced logistic map, cubic map and parametri- 
cally driven Duffing oscillator. In these systems we have 
studied the applicability of 0-1 test to distinguish SNAs 
from quasiperiodic motion and chaos, and have shown 
that the test clearly distinguishes SNAs from other at- 
tractors. The present study has also been able to detect 
the strange nonchaotic dynamics in the Heagy-Hammel, 
fractalization and intermittecy routes, which are in gen- 
eral not easily detected using the Lyapunov exponent. 
We have pointed out that in each of the routes, the 0-1 
test clearly detects an abrupt change in the asymptotic 
growth rate K at critical parameter values during tran- 
sitions. 

In particular for the HH route the value of K is 
smooth in the torus region and at a critical value there 
is a sudden jump after which it exhibits an irregular be- 
havior. In the fractalization route, the value of K grows 
linearly from torus to WT and becomes irregular at the 
transition to SNA. For intermittency route the value of 
K exhibits a sharp transition from a fluctuating value 
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Appendix A: Choosing the parameter c of the 0-1 
test in quasiperiodically driven systems 

In this appendix we illustrate the method by which a 
suitable choice for the parameter c in Eqs.(l)-(2) of the 
0-1 test can be made using the quasiperiodically forced lo- 
gistic map (11) as a specific example. The same analysis 
has been performed for the quasiperiodically driven cubic 
map and Duffing oscillator. Here we present the nature 
of different dynamical attractors in the quasiperiodically 
forced logistic map, which has been studied in Sec. III. B, 
and examine the dynamics in terms of the translation 
variables (p, q) for tori, SNAs and chaotic attractors for 
different values of c chosen in terms of irrational num- 
bers. The results indicate that while many values of c 
fail to distinguish quasiperiodic/SNA/chaotic attractors, 
there exists a set of values which clearly distinguishes 
them as shown in Figs. 13. Particularly, we find that for 

the choice c = ^7-, where ui = (>/5— 1)/2, (see Figs.l3e) 
very clear distinction between different attractors can be 
made. Note that this choice does not produce any res- 
onance with the (p, q) variables. The same value works 
quite well in the cases of quasiperiodically driven cubic 
map and Duffing oscillator. So in our analysis we fix c 

as c = , where u) — (\/5 — l)/2. A second choice 
is c = (V5+ l)/2, corresponding to Fig. 13(a), which we 
find to be reasonably satisfactory, though not as good as 



12 



(iii) 





1 1 
P(n) 


2 




ffpp 



-2.5 -1.5 -0.5 0.5 1.5 2.5 
P(n) 





-250 
-500 
-750 
-1000 



1.5 
0.5 
-0.5 
-1.5 
-2.5 



2 4 6 
P(n) 



-2.5 -1.5 -0.5 0.5 1.5 
p(n) 



1 jf*"*^ 




-1 -0.5 0.5 1 1.5 
P(n) 

800 I ■ ■ — I 

600 \. 
400 

200 \ 

I ^ 



-700 -500 -300 -100 
P(n) 



-4 -2 2 4 
P(n) 



-113 5 7 
P(n) 





-200 
-400 
-600 



10 20 30 40 
P(n) 



-4 -2 2 4 
P(n) 



4 
2 


-2 - 
-2 



-200 
-400 
-600 



2 
P(n) 



10 20 30 40 
P(n) 



10 

-10 
-20 

4 

2 

-2 



-20 -10 10 
P(n) 



-4 -2 2 4 
P(n) 



4 

2 

1 
-2 

-4 
4 

^ 2 
"5" 


-2 



10 
5 


-5 
-10 



2 4 6 
P(n) 



P(n) 



-2 2 6 10 
P(n) 



3 



FIG. 13: (Color online) Choosing the most optimal value of c: Attractors of the quasiperiodically forced logistic map in the 
translation variables (p, q) space, (i) torus, (ii) SNA and (iii) chaotic attractors for the following choices of c: (a) V^> +1 _ 
1.618034 (b) ^" 1 = 0.618034 (c)^( mod 2tt) = 3.883222 (d)^H( mod 2tt) = 5.513002 (e) mod 2tt) = 1.045568 



(f) 2zl( mod 2tt) = 3.883222, where cj = and a = ^4±i Note that the choice (e) with 

most satisfactory description of all the three attractors. Here N — 5 x 10 4 . 

c = 2zeL ! where w = (y/E — 1) /2. 
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Appendix B: Behaviors of the translation variables 
and asymptotic growth rate as a finction of N 



In order to investigate how our results behave as a 
function of the length N of the time series, we discuss 
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dicates torus for e' = 0.3, a = 3.4874; SNA for e = 0.3, 
a = 3.488; chaotic attractor (C) for e = 0.5, a — 3.6. 



for the quasiperiodically forced cubic map and Duffing 
oscillator. 



Appendix C: Characterization of different attractors 
of the quasiperiodically forced cubic map using 0-1 

test 



the dynamics of the translation variables (p, q) as well as 
variation of the asymptotic growth rate if as a function 
of N for the quasiperiodically forced logistic map. We 
have plotted the variables (p,q) in Figs. 14 for different 
types of attractors (Eq.(ll)) for different lengths N of 
time series data. From the figure, one can see first that 
the (p, q) trajectories of the torus motion are bounded 
for N = 5 x 10 4 (Fig.l4.a(i)). Then, upon increasing 
N the (p, q) trajectories continue to be bounded (Figs. 
14a(ii)-14.a(iv)) even upto N — 3x 10 5 . Next, we demon- 
strate in Fig. 15 that the asymptotic growth rate K takes 
a value close to zero irrespective of the length of time se- 
ries (FigT5(a)). We also include the corresponding plot 
of the Lyapunov exponent in Fig. 15(b) for comparision. 
We have also carried out similar analysis for the SNAs 
and chaotic attractors which are reported in Figs. 14(b) 
and 14(c) respectively. The corresponding asymptotic 
growth rate K and Lyapunov exponents are depicted in 
Figs. 15(a) and (b) respectively (Note that in Fig. 15(a), 
the value of K is shown without taking averages over 
segments in order to depict its behavior as a function 
of N). The above results confirm that in the cases of 
SNA and chaotic attractors also the dynamics remains 
qualitatively unchanged irrespective of the length of the 
time series. Finally, the results have also been confirmed 



In this appendix, we present the characterization of 
different dynamical attractors in the quasiperiodically 
forced cubic map (Eq.(12)) in terms of the translation 
variables (p,q) and the mean square displacement M(n) 

using the 0-1 test with c = , and u> = (\/5 — l)/2. 
The results are given in Fig. 1(3. Again this analysis in- 
dicates that the 0-1 test allows one to distinguish SNAs 
from quasiperiodic and chaotic attractors as in the case 
of the quasiperiodically forced logistic map (Sec. III.B). 



Appendix D: Statistical evaluation of the asymptotic 
growth rate and finite time Lyapunov 
exponents(FTLEs) 

In this section we wish to present a statistical evalua- 
tion of our results by analyzing the distribution P(K, N) 
of local asymptotic growth rate of K and comparing the 
same with the distribution P(A, N) of FTLEs for each 
one of the transitions from torus to SNA. The local dis- 
tribution of K is calculated over short segements of the 
trajectories of the given dynamical system. Then, the 
distribution of local asymptotic growth rate of K is de- 
fined similar to FTLEs [21] as 
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FIG. 16: (Color online) For the quasiperiodically forced cubic 
map (Eq.(12)) (a) Torus for / = 0.10, A = 2.14 (b) SNA for 
/ = 0.7, A = 1.86687 and (c) chaotic attractor for / = 0.10, 
A = 2.19: (i) Projection of the attractor; (ii) dynamics of the 
translation variables (p, q) in terms of 0-1 test (Eqs.(l)-(2)); 
(iii) Mean square displacement M(n) as a function of n. Here 
TV = 5 x 10 4 with n = TV/10. 



P(N, K)dK =Probability that K takes a value between 
K and K + dK. 

We have calculated both the distributions of the local 
asymptotic growth rate (P(K,N)) and FTLEs(P(A, TV)) 
in order to compare the effectiveness of the two descrip- 
tions. For each one of the routes to SNAs, we have 
investigated the distributions of local K and FTLES, 
(i) for a torus prior to the transition to SNAs and (ii) 
for a SNA just after the transition. We first present the 
results for the quasiperiodically forced logistic map in 
Fig. 17. The local distribution of K is given in Figs. 17(a), 
17(c), and 17(e) and the distribution of FTLEs is given 
in Figs. 17(b), 17(d) and 17(f) corresponding to HH, 
fractalization, and intermittency routes, respectively. 

From Figures 17(a), (c) and (e), one can note that for 
all the transitions P(K,N) tends strongly to zero when 
the attractor is a torus but on the SNA the distribution 
spreads between 0.1 and 0.9. In particular, the value of 
K is distributed from 0.1 to 0.35 for the birth of SNA 
through HH route. The distribution of K extends from 
0.05 to 0.18 for the logistic map (or 0.2 to 0.3 for the 
cubic map in Fig. 18) for the genesis of SNA through 
fractalization route. In the intermittent SNA the distri- 
bution spreads across the range of K from 0.1 to 0.9. 
The distribution of different ranges of values of K helps 
to distinguish the different transitions to SNAs. It is 
instructive to compare the above results with the distri- 
bution of finite time Lyapunov exponents as shown in 



FIG. 17: A comparative distribution of asymptotic growth 
rate P{K, TV) as a function of K for the three transitions 
showing its values for a torus (A) just before transition and a 
SNA (continuous curve) just after transition in the case of the 
quasiperiodically driven logistic map (Eq. 11): (a) HH route 
for the parameter e = 0.3, a — 3.4874 (for torus) and a = 
3.488 (for SNA); (c) Fractalization route e' = 1.00, a = 2.63 
(for torus) and a = 2.652 (SNA); (e) Intermittency route with 
e = 1.00, and a = 3.40581 (torus), and a = 3.4058088 (SNA). 
The corresponding distribution of P(X, TV) as a function of A 
are given in (b), (d) and (f), respectively, for the above three 
routes for comparision purpose. 

Figs. 17(b), (d),(f). One can note that the distribution of 
finite time Lyapunov exponents have the same shape for 
HH and fractalization route and a slightly different one 
for the intermittent SNA. Moreover, the distribution of 
P(X, TV) for the both cases of torus and SNA are almost 
identical except for the tails on the positive side of A. 
However, the distribution of local K clearly distinguishes 
torus and SNA as well as the three transitions. 

A similar behavior of the distribution of both local 
K and FTLEs has been observed for the cubic map as 
shown in Fig. 18. For the case of quasiperiodically forced 
Duffing oscillator also a similar picture arises. 
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FIG. 18: Same as Fig. 17 but now for the case of quasiperi- 
odically forced cubic map (Eq. 12). The choice of param- 
eters here are as follows. (a),(b): HH route with / = 0.7, 
and A = 1.8865 (torus) and A = 1.8870 (SNA); (c),(d): 
Fractalization route with / = 0.1, and A = 2.14 (torus) 
and A = 2.167 (SNA); (e),(f): Intermittency route with 
/ = 0.7, A = 1.8017 (torus) and A = 1.801685 (SNA). 



[1] P. Grassberger and I. Procaccia, Phys. Rev. Lett 50, 326 
(1993); Physica D 9, 189(1983). 

[2] E. Ott, Chaos in Dynamical Systems, (Cambridge Uni- 
versity Press, Cambridge, England, 1994). 

[3] M. Lakshmanan and S. Rajasekar, Nonlinear Dynam- 
ics Jntegrability Chaos and Patterns, (Springer-verleg, 
Berlin, 2003). 

[4] G. A. Gottwald, and I. Melbourne, Proc. R. Soc. London. 

Ser. A 460, 603 (2004). 
[5] G. A. Gottwald, and I. Melbourne, Physica D 212, 110 

(2005). 

[6] G. A.Gottwald and I. Melbourne, SIAM J. Appl. Dyn. 
Syst 8, (2009); G. A. Gottwald and I. Melbourne, Non- 
linearity 22, 1361 (2009). 

[7] J. Hu, W. Tung, J. Gao, and Y. Cao, Phys. Rev. E 72, 
056207 (2005). 

[8] W. K. Christopher and S. Smith, Phys. Rev. E 83, 03621 
(2011). 

[9] M. Romera-Bastida and A. Y. Rayes-Martinez, Phys. 

Rev. E 83, 016213 (2011). 
[10] C. Grebogi, E. Ott, S. Pelikan, and J. A.Yorke, Physica 

D 13, 261 (1984); C. Grebogi, E. Ott, F. J. Romeiras, 

and J. A.Yorke, Phys. Rev. A 36, 5365 (1987). 
[11] S. P. Kuznetsov, A. S. Pikovsky, and U. Feudel, Phys. 

Rev. E 51, R1629 (1995). 
[12] F. J Romeiras, E. Ott, Phys. Rev. A 35, 4404, (1987). 



[13] M. Ding, C. Grebogi, and E. Ott, Phys. Rev. A 39, 2593 
(1989). 

[14] A. Bondeson, E. Ott, and T. M.Antonsen, Phys. Rev. 

Lett. 55, 2103 (1985). 
[15] J. A. Ketoja and I. Satija, Physica D 109, 70 (1997). 
[16] W. L. Ditto, M. L. Spano, H. T. Savage, S. N. Rauseo, 

J. Heagy, and E. Ott, Phy. Rev. Lett. 65, 533 (1990). 
[17] T. Zhou, F. Moss, and A. Bulsara, Phys. Rev. A 45, 5394 

(1992). 

[18] W. X Ding, H. Deutsch, A. Dinklage, and C. Wilke, Phys. 

Rev. E 55, 3769 (1997). 
[19] G. Ruiz and P. Parmananda, Phys. Lett. A 367, 478 

(2007). 

[20] K. Thamilmaran, D. V. Senthilkumar, A. Venkatesan, 
and M. Lakshmanan, Phys. Rev. E 74, 036205, (2006); 
D. V. Senthilkumar, K. Srinivasan, K. Thamilmaran, and 
M. Lakshmanan, Phys. Rev. E 78, 066211 (2008). 

[21] T. Yang, and K. Bilimgut, Phys. Lett. A 236, 494, 
(1997); Z. Liu and Z. Zhua, Int. J. Bifurcations Chaos 
Appl. Sci. Eng 6, 1383 (1996); 

[22] A. Prasad, V. Mehra, and R. Ramasamy, Phys. Rev. 
E 57, 1156, (1998). A. Prasad, V. Mehra, and R. Ra- 
masamy, Phys. Rev. Lett 79, 4127 (1997). 

[23] A. Venkatesan, M. Lakshmanan, A. Prasad, and R. Ra- 
maswamy, Phys. Rev. E 61, 3641 (2000). 

[24] J. F. Heagy, and H. Hammel, Physica D 70, 140 (1994). 



16 



[25] E. Ott and J. C. Sommerer, Phys. Rev. Lett A 188, 39 
(1994). 

[26] K. Kaneko, Theor. Phys 71, 1112 (1984); Prog. Theor. 

Phys 71, 140 (1994); 
[27] T. Nishikawa, K. Kaneko, Phys. Rev. E 54, 1063 (1996); 

Phys. Rev. E 54, 6114 (1996) 
[28] A. Venkatesan and M. Lakshmanan, Phys. Rev. E 58, 

3008 (1998). 

[29] A. Venkatesan, K. Murali and M. Lakshmanan, Phys. 

Lett A 259, 246 (1999). 
[30] E. J. Ngamga, A. Nandi, R. Ramaswamy, M. C. Romano, 

M. Thiel, and J. Kurths, Phys. Rev. E 75, 036222 (2007); 

E. J. Ngamga, D. V. Senthilkumar, and J. Kurths, Eur. 

Phys. J. Spec. Top 191, 15 (2011); E. J. Ngamga, D. V. 

Senthilkumar, A. Prasad, P. Paramananda, N. Marvan, 

and J. Kurths, Phys. Rev. E 85, 026217 (2012); 
[31] A. Venkatesan and M. Lakshmanan, Phys. Rev. E 63, 

026219 (2001). 

[32] U. Feudel, S. Kuznetsov, and A. Pikovsky, Strange 
Nonchaotic Attractors: Dynamics between Order and 
Chaos in Quasiperiodically Forced Systems, (Worls Sci- 
entific, Singapore, 2006). 

[33] A. Prasad, S .S Negi, and R. Ramaswamy, Int. J. Bifurc. 
Chaos. Appl. Sci. Eng 11, 291 (2001). 

[34] S. Y. Kim, W. Lim, and E. Ott, Phys. Rev. E 63, 056203 
(2003). 



[35] Y. C. Lai, Phys. Rev. E 53, 57 (1996). 
[36] T. Yalcinkaya and Y. C. Lai, Phys. Rev. Lett 77, 5039 
(1996). 

[37] X. Wang, M. Zhan, C. H. Lai, and Y. C. Lai, Phys. Rev. 

Lett 92, 074102, (2004); X. Wang, Y.C Lai, and C. H 

Lai, Phys. Rev. E 74, 016203 (2006). 
[38] A. Khavanov, N. A Khavanova, P. V. E Mcchintock, and 

S. Anischenko, Phys. Lett. A 268, (2000). 
[39] U. Feudel, J. Kurths, and A. Pikovsky, Physica D 88, 

176, (1995); U. Feudel, C. Greboji, and E. Ott, Phys. 

Rep 290, 11,(1997). 
[40] A. Pikovsky and U. Feudel, Chaos 5, 253 (1995). 
[41] S .S. Negi, A. Prasad, and R. Ramaswamy Physica D 

145, 1 (2000). 

[42] M. Nicol, I. Melbourne, and P. Ashwin, Nonlinearity 14, 
275 (2001). 

[43] P. Ashwin, I. Melbourne, and M. Nicol, Physica D 14 
275 (2001). 

[44] J. H. P. Dawes and M. C. Freeland, The 0-1 test 
for chaos and strange nonchaotic attractors; peo- 
ple. bath. ac.uk/jhpd20/publications (2008); 

[45] W. H. Press, S. A. Teukolsky, W. T. Vetterling, and 
B. P. Flannery, Numerical Recipes in C (Cambridge Uni- 
versity Press, Cambridge, 2002). 



